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Chapter 1 


Introduction 


The great charm resulting from this consideration (a multi-connected universe) 
lies in recognition of the fact that the universe of these beings is finite and yet has 
no limits. 

- Albert Einstein, on a multi-connected universe |1; 2] 


In the year 1916, Einstein hypothesized the general theory of relativity, which 
exposed gravity as a manifestation of the curvature of space-time and further 
relegated the extent of curvature to the distribution of the surrounding mass. 
However, although this theory made it possible to predict the local geometry of 
space, it did not put any constraints on its boundary conditions or the global 
geometry of space-time. 


Even though we encounter the significance of the spatial curvature in cosmol- 
ogy texts relatively often, the possibility of a multi-connected topology is usually 
not discussed. However, that does not imply that the possibility does not exist. 
Einstein’s first seminal paper in cosmology discussed the possibility of a multi- 
connected universe according to a Clifford-Klein space form [3]. Einstein and 
Wheeler favored a finite universe on the basis of Mach's principle [4]; others have 
said that an infinite universe is unaesthetic and wasteful [5], while quantum cos- 
mologists say that small universes have small action and are thus more likely to 
be created [6]. 


On different distance scales, the geometric representation of space changes, 
for instance, for distances between 10718 and 101! metres, it can be approximated 
by ordinary three dimensional euclidean space, E?, while on larger scales of up to 
about 107° metres, space is better described by a continuous Riemannian mani- 
fold where the curvature varies in response to the surrounding mass content. 


Topology is essentially a study of continuity of surfaces and manifolds. It 
is a property which is preserved through kneading, stretching or squeezing the 
manifold, actions which change the metric. However, cutting, tearing or making 
handles in the manifold would result in a different topology. A famous example 
would be of a doughnut and a coffee cup, which are equivalent entities accord- 
ing to a topologist, since they can be continuously transformed into each other. 
These kind of transformations are called homeomorphisms and are discussed in 
detail in Chapter 4, along with other topological entities useful for classifying the 
universal manifold. 


The first set of constraints that can be applied to limit the large number of 
possible topologies are homogeneity and isotropy, as is confirmed through obser- 
vations which point at statistical isotropy and homogeneity in the distribution of 
astronomical objects [7]. However, a special class of models also exist that sus- 
pend the assumption of isotropy, Bianchi universes, which I discuss in Chapter 
8. A better understanding of the best possible spaces that model our universe 
can be obtained by studying the vast range of abstract spaces available to us 
from topology. Cosmic Topology thereby combines these two areas, topology and 
cosmological observations, to form a better picture of what the universe looks like 
on scales comparable to its size. 


The most promising source of cosmological data for the purpose of topology 
at the moment is the Cosmic Microwave Background, an ancient relic of the Big 
Bang [8]. With the advent of improved technology, we have been able to send 
satellites specifically for the purpose of measuring the microwave background, pro- 
viding us crucial information on the content and structure of the universe about 
300,000 years post big bang. This radiation contains tiny temperature fluctu- 
ations that reveal the pattern of density fluctuations in photon gluon plasma, 
before the photons decoupled and went on to stream freely through space, form- 
ing the CMB. These fluctuations were first detected by the Cosmic Background 
Explorer (COBE) in 1992 [9], for which George Smoot and John Mather received 
the 2006 Nobel Prize in Physics. High precision measurements were further made 
with the launch of another probe, the Wilkinson Anisotropy Probe (WMAP) in 
2001 [10], which was 45 times more sensitive and had a 33 times higher resolution 
than the COBE satellite. The measurements were recently further improved with 
the Planck satellite [11], which was 10 times more sensitive and had 3 times the 
resolution of WMAP and mapped the sky in 9 instead of the previously used 5 
frequencies. 


The results obtained from the probes largely confirmed the predictions from 
the ACDM model (flat universe model) of the universe except for some anomalies 


in the large scale measurements. The quadrupole and the octopole were mea- 
sured to be much lesser than the value predicted by the model and this has been 
linked to a possible compactification of dimensions in one or two directions in 
space [39]. This can be understood in analogy to a tied string which allows only 
a certain number of wave-modes to exist with a limit on the highest mode being 
the length of the string. As different harmonics come together to form a musical 
note, in a similar way the CMB is a combination of spherical harmonics of tem- 
perature fluctuations. The modes that are sustainable in an instrument give a 
lot of information about the shape of the instrument. For instance, weak second 
and fourth harmonics are characteristic of a clarinet. Hence, for a continuously 
defined field in space, we can split it into a sum of harmonics. These harmonics 
are eigen-modes of the Laplace operator, and in principle, vibrational modes of 
space. Hence just as the density fluctuations in a clarinet are a byproduct and 
hence a characteristic of its size and shape, the pressure and density fluctuations 
in the primordial plasma were constrained with the size and shape of the uni- 
verse. Though we do not see the complete 3D modes we see their intersection 
with the 2D horizon sphere, and it is not hard to get the 3D modes from their 
2D manifestations. 


Some of the suggested manifolds currently in focus, which incorporate some of 
the observed anomalies, are the Poincare Dodecahedron, the Picard horn and the 
Bianchi universes. We also discuss three dimensional manifolds of constant cur- 
vature and the observational consequences of Euclidean, Spherical or Hyperbolic 
manifolds. A central premise of the thesis is the classification and properties of 
different topological manifolds and their relevance to our universe. 


Chapter 2 


Standard Cosmology 


Cosmology essentially focusses on the study of the origin, structure and evolu- 
tion of the universe. A central assumption in cosmology is that the universe is 
isotropic and homogenous at large scales, hereby referred to as the ‘the cosmologi- 
cal principle’. It is however difficult to observe at what scale the premise becomes 
absolutely true, since even at large scales we observe small deviations from ho- 
mogeneity and isotropy which are known to be seeds from which the large scale 
structures later emerged. The seeds themselves are believed to have originated 
from initial quantum fluctuations which were greatly magnified during the era of 
inflation. However, it would be safe to assume that the universe is statistically 
homogeneous at scales larger than about 100 Mpc. 


Relativistic cosmology gives us an idea of the local properties that a manifold 
must possess as deduced from the Einstein’s field equations. However, these equa- 
tions, being partial differential equations, do not say anything about the global 
properties of space-time. Hence, a given local metric element could correspond 
to a large variety of topologically distinct universe models. Below are described 
some core concepts of cosmology which will be later utilized in subsequent chap- 
ters. More detailed discussions can be looked up in Peacock’s [12] and Kolb’s [13] 
classic texts on cosmology. 


2.1 The Friedmann-Lemaitre-Robertson- Walker 
Model 


A homogeneous and isotropic universe is one that can be sliced into maximally 
symmetric 3 spaces of constant curvature and these symmetries constrain greatly 
the possible allowed solutions to the global metric. Friedmann, Robertson and 


Walker proposed a metric which gives us all possible solutions for such constant 
curvature universes. The solutions include big bang solutions, de Sitter solutions 
and also includes those requiring a cosmological constant. 


The Friedmann-Lemaitre-Robertson-Walker metric is of the form: 


dr? 


ds? = g,,dx"dx" = —dt? + a?(n)| + r^ (d? + sin?6d¢’)], (2.1) 


1— kr? 


where a(t) is the scale factor and & — 0, —1,1 corresponds to a flat, positively 
curved or negatively curved universe respectively. 


Converting to conformal coordinates where the manifold can be thought of as 
a constant curvature manifold and the dynamics is incorporated into the confor- 
mal scale factor a(7) using the definition of conformal time dr = dt/a(t). The 
expression now becomes 


ds? = a? (n) |- dr? + do’), (2.2) 
where the spatial part of the metric is given by 
do? = dx? + f (x)(d6? + sin?6d¢’), (2,3) 
f(x) is a function of curvature given by 
y? r=x flat, 


f(x) =< sinh?xy r=sinhy hyperbolic, (2.4) 
sin?y r=siny spherical. 


Hence the FLRW metric of space-time can be expressed as 


—1 0 0 0 
[j^ == 0 0 

jw — | o ur ape 0 (2.5) 
0 0 0  dr?sin?0 


Calculating the Einstein tensor G$ from the metric gives 


Gs = =o 8 m 35) (2.6) 
E :2 
0 u a a K 
pues -(22 +54 ss 2.7 


Substituting this along with the energy momentum tensor T”, into the Ein- 
stein’s equations will give us a system of equations relating different cosmological 
parameters. The energy momentum tensor for a maximally symmetric space-time 
has a simple form given by 


(2.9) 


COO?) o 
Ors Oo 
"S OO. © 


where p denotes the energy density and p the pressure and these quantities only 
depend on the time in a homogeneous universe. 


2.2  Einstein's General theory of Relativity 


Einstein related the curvature of space-time to matter with the local Einstein 
equations: 
Guy = 81GT,,, (2.10) 


where G',,, represents a function of the metric of spacetime and T,,, is the energy 
momentum tensor representing the matter distribution in space. These equations 
quantify the influence of matter fields on local curvature, however they do not 
determine the large-scale topology of the universe. 


The Einstein's equation G,,, = 81GyT,, is a non-linear system of ten partial 
differential equations. In the case of a Friedmann-Lemaitre-Robertson-Walker 
(FLRW) universe, it reduces to two ordinary differential equations which can be 
rearranged to give us the Friedmann equations below: 


Sgt 93 = 8mGwp, (2.11) 
32 = -A4nGw(p- 3p). (2.12) 
a 


From these equations, one can obtain the energy continuity equation given by 


p=-3(p +p)" (2.13) 


a 

From Equation 2.11 we deduce that the global energy density is a function 
of curvature. The correspondence is made in terms of the density parameter, 
Q = p/p. where p, is the total density of a flat universe and can be obtained by 
taking « = 0 in eqn. 2.11, which gives us p, = 3H?/8nG. 


Writing 2.11 in terms of Q, we get the expression H?a?(Q — 1) = &, where it 
becomes evident that Q > 1 for k = 1 and Q < 1 for k = —1. Hence the universe 
will be positively curved for k = 1 and negatively curved for k = —1 [14] 


In terms of the curvature radius, (Reury = a/|k|/?), we have 


1 


Pico = TES a4nNH/9* 
Ha -1[72 


(2.14) 


2.3 The Last Scattering Surface 


About 300,000 years after the big-bang, the universe was a hot plasma of coupled 
baryons and photons. As the plasma cooled down to below the photon tem- 
perature, the photons decoupled from the mixture and started streaming freely 
through space thereby making it transparent. Thus, the surface of last scatter 
is the farthest that we can see into the universe, since that is the time when 
the photons decoupled from matter and were free to stream through space. So 
for a topology to be observable with current methods, a manifold's fundamental 
domain would have to be smaller than the diameter of the surface of last scatter 


(SLS). 


It would thus be useful to calculate the radius of the surface of last scatter, 
which is the distance light has travelled from the time of decoupling to today, 


D, = a(t) T L3 (2.15) 


or in comoving units, 


The diameter of the SLS, 2A7 essentially represents the extent of the observ- 
able universe. 
The volume of the observable universe can thus be determined by integrating the 
metric of space-time in radial coordinates, 


Vsrg = Jv —gdrd0do. 


For Rgrs = An in FLRW space-time metric as given in Equation 2.4, The 
corresponding volume of the observable universe is 


Vsrs = m(sinh(2Am)—2Am) OQ <1, 
4 f 
= An’, Qe 
= m(2Am-— sin(2Am), 909» 1. (2.16) 


For multi-connected spaces, the number of copies of the fundamental domain 
that can fit inside the observable universe would be given by [14] : 


Vsrs 
N= . 
Vu 


where Vsrs is the surface of last scatter and Vm is the volume of the fundamental 
domain. 


2.4 The Cosmic Microwave Background 


The Cosmic microwave background though having being predicted since the 
1940’s was only accidentally discovered in 1964 by the American radio astronomers 
Arno Penzias and Robert Wilson who were later awarded the Nobel prize for the 
discovery. Probes were later sent to map the microwave background, namely 
COBE, WMAP and recently the Planck Satellite. The latest high resolution 
map of the CMB was release by the Planck satellite on 21st March, 2013 as 
shown in Figure 3.1 below. 


The CMB is largely isotropic with tiny fluctuations which are now known to 
represent density perturbations in the early universe. The photons that produce 
the tiny fluctuations in the temperature spectrum are supposed to have origi- 
nated at the surface of last scatter where they decoupled from baryons. The 


Figure 2.1: CMB temperature anisotropy map from Planck satellite. 


temperature fluctuations are made up of two components: the original temper- 
ature fluctuations in the baryon-photon fluid at the time of decoupling and the 
additional component introduced while on its journey to where we are today, 


ôT ôT ôT 
T) obs i T) intr i T). n 


The Temperature perturbations ôT can be defined either as the deviation 
from the averaged out temperature of the background or just the background 
temperature. The latter case is simpler and doesn't take the contribution of the 
fluctuations into the mean, though it is more common to use the former method. 


'The CMB radiation is uniform across the sky and exhibits small temperature 
fluctuations of 1 part in 10°. The temperature fluctuations can find their source in 
the frozen acoustic waves in the baryon photon plasma as the photons decoupled. 
The fluctuations in the photon temperatures arose from escaping from different 
density regions in the plasma, the photons from the denser regions lost more en- 
ergy in escaping the plasma than the ones from rarer regions. The photons, in 
their journey towards us, became redshifted as well. 


The temperature anisotropy can be expressed as a function on a sphere. In 


analogy to the fourier expansion in three-dimensional flat space, we can express 
the fluctuations as a expansion of spherical harmonics as 


ôT - 
l=1,m=-l 


The multipole coefficients, aj, can be calculated from: 
óT 
aw = [ Yi,(0.0)77.0.0)40. (2.19) 


We take ago = 0, since Yoo is a constant and we are left with the average 
temperature. The l = 1 part represents the dipoles due to the motion of the solar 
system with respect to the last scattering surface interspersed with the actual 
cosmological dipole arising from large scale perturbations. Since it hasn't been 
possible to separate the two dipoles, the analysis of the CMB is taken from / = 2 
onwards. 


The multipole coefficients represent a deviation from the average temperature 
and are gaussian in nature having an expectation value of zero. The variance of 
the coefficients, (|ajm|*) averaged over m gives us the average deviation. True to 
the isotropic nature of the background, we only observe dependence on 1, which 
corresponds to the angular size of the mode and none on m, which represents 
orientation. Thus we have 


1 
Q = —— pele 2.20 
= py Lila (2.20) 
where Cj is the angular power spectrum and contains all statistical information 
about the CMB. We compare the power spectrum predicted by different metric 


solutions to the one observed as a test of its validity. 
The observed Cj; is the manifestation of a specific instance from an array of 


random possible outcomes. The observed angular power spectrum is the average 
of the observed am values, 


10 


ao 1 2 
Ce pP lai |?. (2.21) 


m 


On larger scales, the CMB begins to be dominated by the Sachs-Wolfe effect. 
It is given by 


ôT 2 o 
(F) = —=O(taec, Lis) + Diis Tis) + 2 | dt (2.22) 
T obs 3 dec 
1 So: 
= 39 aee, Tis) + a dt. (2.23) 
dec 


The first term in the above equation, SO (Lacs, Lis) is the ordinary Sachs Wolfe 
effect and is thought to have originated due to a gravitational redshift at the 
surface of last scattering. The integrated Sachs-Wolfe effect on the other hand, 
2 ie $dt originated between the surface of last scatter and the earth due to the 
presence of energy densities other than matter in the universe, for instance ra- 
diation energy or dark energy and thus is negligible in a universe dominated by 
matter density. 
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Chapter 3 


Basic Conditions required of a 
universal manifold 


We get a set of solutions from the general relativity equations which represent 
possible cosmological models. We can represent these models as an ordered pair 
(M, g) where M is a connected smooth, n-dimensional manifold (n > 2) and 
g is a smooth Lorentzian metric on M. The information on the topology of the 
universe is contained in M while g relates different points on the Manifold, M. 
The metric, g can be thought of as the function which assigns length to each 
vector £^ in M, classifying them into three classes: timelike if g,,€¢€° > 0, light 
like if g,,£^£^ = 0 and space like if g,,£"£^ < 0 and gives them positive, zero and 
negative lengths respectively. The lightlike curves form a double cone structure 
at every point on the manifold, encompassing all timelike curves. The manifold, 
must admit a timelike vector field for a lorentzian metric to exist. A detailed 
analysis of causality conditions and further lower level constraints on a space- 
time manifold are discussed in [15; 16] 


3.1 Time Orientable Manifolds 


The world lines for a particle are represented with a cone like structure with 
two lobes for past and future oriented vectors. Special Relativity holds that a 
particle can only move on its worldline into the future and thus preserves the 
local orientability of time. This property holds also in curved spacetimes where 
general relativity takes over, with a straightforward extension of the local validity 
of special relativity. However, for a consistent global time orientation definition, 
local time orientations must be continuous and not clash, especially so, in closed 
trajectories, where time orientation must remain the same after a complete trip. 
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Figure 3.1: A double cone structure on a manifold with timelike, null and space 
like vector in, on and outside the cone respectively. 


A simple test to check time orientability of a curve is to consider a point P 
lying on a curve C with an orientation along one of the direction, now take the 
point along the curve such that it returns to its previous position. If the orienta- 
tion is still preserved, the curve is time orientable. A manifold is said to be time 
orientable, if every closed curve on the surface is time orientable. 


3.2 Causal Manifolds 


Notion of causality derives from the fact that cause precedes effect. We say that 
a point p € M chronologically precedes q € M, i.e. p « q if there is a timelike 
path from p to q. We can then define the future cone of p as 


It(p) ={qEeM:p <q}. (3.1) 


For some interval, I C IR, a smooth curve y : / > M is timelike if its tangent 
vector £^ at each point in y is timelike. Similarly, a curve is null if its tangent 
vector at every point is null. A curve thus is causal if its tangent vector at each 
point is either null or timelike. For a space-time manifold to be physically pos- 
sible, it should not have any closed timelike curves i.e. points which violate the 
chronology condition listed above. A manifold is further called strongly causal if 


13 


for each p € M there exist arbitrary small neighbourhoods U of p such that any 
causal curve that starts in U and then leaves U never re-enters U. 


In special relativity, it implies that a real particle cannot travel along closed 
spacelike curves, which is another term for travel into the past. However in gen- 
eral relativity, it may be possible to have a view of the past while remaining in the 
future light cone, due to subtle distortions of spacetime by strong gravitational 
fields from a rotating massive body. 


To find a stable solution, we can define stably causal manifolds. A spacetime 
is stably causal if it admits a cosmic time function that is a continuous real func- 
tion and whose gradient is universally timelike. Some stably causal spacetimes are 
Minkowski, Schwartzchild, Friedmann, which are globally time orientable mani- 
folds where time must increase along causal worldlines. 


3.9 Global Hyperbolicity 


A Cauchy surface S for M is an achronal (i.e. no timelike paths exist between any 
two points) hypersurface S in M which is met by every inextensible causal curve 
in M. Or simply put, it is the hypersurface constituting all of space at a given 
instant of time. Further, if S is a Cauchy surface for M, then M can be shown to 
be homeomorphic to IR x S. 

So in the simplest terms, Globally hyperbolicity can be defined as being equiva- 
lent to the existence of a Cauchy surface or the ability to foliate the manifold with 
Cauchy surfaces. Such a manifold is necessarily stably causal and time oriented. 


The concept of determined path of evolution of a system from a certain point in 
space can extended to spacetime through the concept of domain of dependence. 
For an initial spatial hypersurface X, the domain of dependence is a region of 
spacetime D*(X), such that any timelike curve reaching any point in D*(*) 
must intersect €. Similarly D^ (X) is the past analog of D*(=). Hence we have 
a region of space-time 

DOS Sb" p Os), 


which is entirely determined by information on X. Imposing this condition on 
spacetime infers that it is possible to determine the structure of spacetime from 
information on a single hypersurface, analogous to an initial value problem. It 
also implies that the interior of M can be taken to be diffeomorphic to IR x Ms 
where M is a 3 dimensional Riemannian manifold. 

A manifold without boundary is said to be globally hyperbolic if two conditions 
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hold: 

1. For every pair of points p and q in the manifold, the space of all points that 
can be reached from p along a past oriented causal curve and from q along a 
future oriented causal curve, is compact and can be denoted by 


D^ (p)  D* (q). 


2. Causality holds on the manifold. 


3.4 Spatial Orientability 


Additional restrictions could be imposed on spacetime based on space orientabil- 
ity. Space orientability can be easily described in terms of a two dimensional 
surface. On taking an arrow around the surface, if it returns to its initial position 
and still points in the same direction, the space is orientable and if not then it is 
non-orientable. A well known example of a non-orientable surface is the Mobius 
strip, which reverses the direction of the normal to the surface on complete trip 
around the strip, see Figure 3.2. 


Generalizing to higher dimensional manifolds, we can define two types of 
curves, left handed and right handed. In three dimensions, a person crossing a 


non-orientable manifold, would have their left and right sides switched. Thus in 
general terms, a manifold which preserves the handedness of a closed curve is 


spatially orientable. 
n 


Figure 3.2: A Mobius strip is a good example of non-orientable 2-dimensional 
manifold. [15] 
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Chapter 4 


Topological Classification of 
Spaces 


Topology can be defined as the study of continuous transformations, albeit the 
property which remains unchanged when continuous transformations are made 
on a geometry, like squeezing, stretching etc. which change its metric but not 
its topology. Two manifolds that belong to the same topological class are called 
homeomorphic and can be continuously and reversibly transformed into each 
other. In other words, if we have two manifolds Mj and Mə, a homemorphism 
between them would be continuos map ® : M; — M» with a defined inverse. 


FLRW Models admit spatial sections of homogeneous and isotropic spheri- 
cal, hyperbolic or euclidean manifolds depending on whether the sign of spatial 
curvature is positive, negative or zero. However, there is often a common mis- 
conception that flat or hyperbolic universes imply an infinite universe, which was 
proved unfounded long ago by Friedmann [17] and Lemaitre [18] who discovered 
that FL metrics with zero or hyperbolic topologies did admit spatially closed 
topologies. However, the works of these and many other people remain ignored 
and cosmology textbooks implicitly assume space to be a simply connected hy- 
persphere. 


Given the isotropy of the microwave background, it is implied that the curva- 
ture of space is almost constant throughout. Hence literature on possible mani- 
folds for the universe focuses mainly on manifolds of constant curvature. General 
relativity is invariant under diffeomorphisms which signify change of coordinates 
but not homeomorphisms. Thus, the principle of covariance goes a long way in 
predicting the laws of physics that a body follows in different regions of space, 
however topology remains independent of such a correlation. 
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Figure 4.1: Some examples of homeomorphic manifolds, the numbers below the 
columns refer to the number of holes in the topology [15]. 


Below are explained some basic concepts that are frequently used in classifying 
and studying topological spaces: 


4.1 Homogeneous Spaces 


For two dimensional surfaces, it was shown by [19], that if a space is closed and 
connected, it is homeomorphic to Riemannian surfaces of constant curvature. A 
Riemann surface is defined as a complex manifold of dimension one. Hence all 
closed surfaces can be classified into one of the three Riemannian metrics: Spher- 
ical S?, Euclidean E? and Hyperbolic H?. 


In three dimensions, this is not true, as we can clearly see from a three di- 
mensional cylinder S? x R which is not homeomorphic to any of the constant 
curvature geometries, neither the spherical nor the euclidean manifold. The met- 
ric of the three dimensional cylinder is homogeneous but anisotropic. There are 
in total eight types of homogeneous three dimensional geometries out of which 
only five of them are of a constant curvature. 


The symmetries of a manifold can be quantified with the group G of isometries, 
which are transformations to the manifold that leave the metric invariant. For 
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a homogeneous manifold, G is non-trivial. The group H of all points y acts 
transitively on M as it relates to g € G such that g(x) = y with y being referred 
to as the orbit of z. The subgroup of isometries that leave the point x fixed (for 
instance, rotations around x) is called the isotropy group, I at x. These isometries 
are related by 

dim(G) = dim(H) + dim(1). (4.1) 


G is simply transitive on H if dim(G) = dim(H) and multiply transitive if 
dim(G) > dim(H). 


The isometry group has dimensions < n(n + 1)/2 for an n-dimensional man- 
ifold and attains maximum value for a maximally symmetric manifold. A maxi- 
mally symmetric manifold is essentially a manifold which has the same number 
of symmetries as an ordinary euclidean space. For the space-time metric of our 
universe, a maximally symmetric space should have an isometry group of dimen- 
sion = 10. 


As [15] analyses, space times of dimensions (> 6 < 10) are not realistic cos- 
mological models due the large number of dimensions involved. For dim(G) < 6, 
the group may act on M or lower dimensional manifolds. For dim(G) = 4, the 
isometry acts on a manifold that is homogeneous in space and time and gives us a 
model which is spatially symmetric in space and time, but does not allow expan- 
sion. For the alternative case, when a subgroup of G acts only on the space-like 
hypersurfaces, giving us spatially homogeneous space times. This scenario has 
three possible sub cases: 


- dim(G) = 6, where G3 acts on spatially homogeneous spaces and there is a 
G3 isotropy group. Thus, we have spatially homogeneous and isotropic space 
times which also allow space like hypersurfaces of constant curvature i.e. the 
FLRW models. 


- dim(G) = 3, where there is only one group of isometries, i.e. the group Gs 
acting on the spatially homogeneous spaces. Thus we have homogeneous and 
anisotropic spaces, one example of which are the Bianchi models [20]. 


- dim(G) — 4, where G is multiply transitive on 3 dimensional subspaces, some 


such space times are discussed in [21]. We will not consider these space-times in 
this text. 
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4.2 Simply and Multiply Connected Spaces 


A good place to begin is the definition of the concept of homotopy which is an 
important classification regime used in topology. Two loops y and ^" drawn on 
a manifold surface are said to be homotopic if one can be continuously trans- 
formed to the other. A simply connected manifold now can simply be defined as 
a manifold for which any loop is homotopic to another, or equivalently, all loops 
are homotopic to a point. If this condition is not true for all possible loops on 
the manifold, it is multi-connected. Homotopic loops give us information about 
holes or handles in a manifold. In higher dimensions, one dimensional homotopy 
loops are not enough to encompass all the properties of the topology, leading to 
the introduction of homotopy groups. The first homotopy group is called the fun- 
damental group. Poincare conjectured that any connected closed n-dimensional 
manifold with a trivial fundamental group is topologically equivalent to a sphere 
[22]. 


Multiconnectedness implies that the fundamental group is non-trivial, essen- 
tially meaning that there is one hole in the manifold. Poincaré in 1904, conjec- 
tured that a connected closed n-dimensional manifold with a trivial fundamental 
group must be topologically equivalent to a sphere, S”. 


The set of solutions to Einstein's equations does not place any topological 
constraints on the manifold except its curvature. The FLRW models describe the 
observed universe with the greatest accuracy among the known models and give 
solutions for homogeneous and isotropic models with spherical, hyperbolic or flat 
topologies further incorporating a wide variety of possible solutions like the de 
Sitter solution or solutions with a cosmological constant or a non standard equa- 
tion of state. The assumption in most literature of a simply connected universe 
is arbitrary and replacing the same with a multi-connected universe changes a 
very few characteristics in the FLRW models. One of the differences lies in the 
range of the coordinates where for a simply connected universe, one would have 
$:0— 27,0: —1/2 > 1/2, x :0— œ for k = —1,0 and 0 —^ m fork = 1 
whereas for a multi-connected model, space is smaller and the range of the coor- 
dinates is reduced. 


As discussed earlier, observations suggest that the universe is homogeneous 
and locally isotropic which implies that space has constant curvature. Thereby, 
most multi-connected models explored in literature rely on this assumption with 
the exception of Bianchi and Lemaitre-Bondi models among a few others. How- 
ever, even with the anisotropic models, the homogeneity and local isotropy of 
these models ensure that the CMB remains isotropic. A significant difference 
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however is observed in the spectrum of density fluctuations. 


While the finiteness of simply connected models can simply be determined 
from the sign of curvature of the manifold, i.e. infinite for k = 0, —1 and finite for 
k — 1, the same does not hold true for multi-connected topologies. As early as 
1924, it was known that multi-connected models with a zero or negative curvature 
admitted spatially closed topologies [23; 24]. For instance a toroidal universe is 
of a finite volume and circumference despite zero curvature. 


4.3 Fundamental Domain 


A simple example is a torus whose fundamental domain is a rectangle. To obtain 
a torus from a rectangle, we first identify one pair of opposite sides in the rect- 
angle, thereby getting a cylindrical tube. Identifying the other pair of opposite 
sides gives us a torus. 


The transformations done in identifying the opposite edges form a holonomy 
group. A holonomy group is a subset of the full isometry group of the covering 
space. To understand the holonomy group, consider a point x and a loop y at 
x in the manifold M. If y lies in a simply connected domain of M, it generates 
a single point Z in M but if the manifold is multi-connected, it creates a set of 
points 2", Z"... which are said to be homologous to z. The displacements form 
isometries referred to as the holonomy group I in M. Since there is a non zero 
distance between the homologous points, the group is discontinuous and has no 
fixed generating point. Thus, the holonomy group is said to act freely and dis- 
continuously on M [15]. An example of homologous points on a 2 torus manifold 
is shown in Figure 4.2 


The full isometries of the universal cover are broken by identifications and we 


can represent a compact manifold as a quotient space given by G/T where G is 
the group of isometries of the domain and I is the holonomy group. 


4.4 Universal Covering Space 
The universal cover of a connected topological space is a simply connected space 


with a map £ Y — X that is a covering map. By acting with the transfor- 
mation group on the fundamental domain, we get many identical copies of the 
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Figure 4.2: A representation of holonomic generators in a 2-torus manifold. 


domain which give us the universal covering manifold. In case of simply con- 
nected spaces this universal covering is identical to the fundamental domain for 
instance a sphere S? is its own universal cover, however in the case of multiply 
connected spaces we get replicas of the central manifold. In this case, a univer- 
sal manifold is constructed as follows. The space is cut open to make a simply 
connected space with edges, called the fundamental domain of the manifold. For 
instance, a hyperbolic octagon for a double torus and a square for a square torus. 
Now add another copy of the fundamental domain to the edge and keep doing so 
until all edges of the original manifold are covered. More copies are added to the 
resulting space recursively until a covering map with possibly infinite number of 
copies of the fundamental domain is obtained. The largest such possible cover is 
called the universal covering space. Thus, If f : Y — X is a covering map, then 
there exists a covering map f: X— Y such that the composition of X and X is 
the projection from the universal cover to X [25]. 


For example, a flat torus tiling a universal covering space can be likened to 
the screen of video game where on walking off the right end of the screen, one 
would emerge on the opposite edge and the same for the vertical edges. Thus, one 
gets the impression of an infinite space even though it is just a repetition of the 
same fundamental domain over and over. This type of a universal covering space 
is constructed by identifying the edges of a fundamental domain, and identifying 
the edges differently gives rise to a different set of orientations and symmetries 
on transition between different universes. An example of a universal covering of 
a 2-torus is given in Figure 4.4 
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Figure 4.3: The universal covering space of a 2-torus, the space is tiled with 
squares which bear specific orientation and symmetry relations to each other 
[14]. 


4.5  Detectability of a Multi Connected Topol- 
ogy 


There are basically three possible correlations between detectability of topology 
and the size of the universe. First, that the universe could be infinite in which 
case, it is not possible to detect topology with currently known methods. Sec- 
ondly the universe is finite but much larger than the scale of the observable 
universe in which case, it is hard to detect visible signs of topology. The last and 
the best scenario would be a universe which is finite and comparable to the size 
of the observable universe, where we can thus use current methods discussed in 
Chapter 8 to detect its structure. 


For a manifold M, there can be defined the smallest and largest circles in- 
scribable in M, Tinj, described in terms of the smallest closed geodesic, lm and 
Tmax respectively. A closed geodesic, that passes through a point x, in a multiply 
connected manifold is a part of the geodesic that connect that point to its image 
in the covering space M. The length of any such closed geodesic which passes 
through x, in a manifold with a fixed isometry g, is given by the distance function: 


óg(x) = d(x, ux). (4.2) 
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In terms of this distance, the injectivity radius can be defined as 


rwr) = Zminger{ ögla), (4.3) 


where I' denotes the covering group without the identity map. 


We can define the observation survey depths to be Xobs. A topology is said 
to be undetectable if Xobs < Tinj in which case we cannot detect any multiple 
images in the observable sky, and detectable for Xobs > Tinj. 


Figure 4.4: Two fundamental cells and the smallest and largest sphere's inscrib- 
able in M [26]. 


In a globally homogeneous manifold, the distance function for any covering 
isometry is constant, hence fin; is constant throughout space and the detectability 
of a topology does not depend on the observers position in the manifold, In 
inhomogeneous manifolds however, r;,; varies from point to point and thus the 
topology depends on both the observers position and the survey depth. In this 
case however, we can still define an absolute undetectability condition, that is if 
for 

Tinj = Mingem{ Tinj(£)}, 


Xobs < Tinj 
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then the topology is undetectable for all observers in M [26; 27]. 


For a flat manifold E?, the relationship between the horizon radius and the 
injectivity radius is pretty arbitrary since it is possible to stretch its translational 
components to obtain any injectivity radius. It is thus probably least likely to be 
able to detect a Euclidean topology from the three possible topologies. On the 
other hand, in the case of a Hyperbolic topology, the volume of the fundamental 
domain increases as the complexity of the Hyperbolic group I increases. However, 
the injectivity radius for the smallest hyperbolic manifolds exceeds the horizon 
radius. Spherical manifolds on the other hand decrease in size as the symmetry 
group I becomes larger. 


24 


Chapter 5 


Three Dimensional Manifolds of 
Constant Curvature 


Current observations of the observable part of the universe imply a homogeneous 
and isotropic geometry to a precision of 1 part in 104, thus we should begin by 
considering topologies that are locally homogenous and isotropic, thereby the 
constant curvature Euclidean, Spherical and Hyperbolic geometries as deduced 
in the previous chapter. 


Any compact 3-manifold M with a constant curvature k allows a discrete 
isometry subgroup I acting freely and discontinuously on M. Such a manifold 
can also be written as M/T, where M is the universal covering space of M, given 
by Euclidean space (for k = 0), 3-sphere (for k > 0) or the hyperbolic 3-space 
(for k « 0). [15] does an excellent work of classifying the Euclidean, Spherical 
and Hyperbolic manifolds further into sub manifolds as discussed below. 


5.1 Euclidean Manifolds 
The line element of the Euclidean covering space is given by: 
do? = R?( dx? + X?(d6? + sin?0dq?))). (5.1) 


where x > 0. 


The full isometry group is given by G = ISO(3) = R? x SO(3) and the genera- 
tors of the possible holonomy groups I include different combinations of identity, 
translations, glide reflections and helicoidal motions. In total 18 different types 
Euclidean manifolds can be generated. The manifolds can be classified primarily 
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into open and closed models. 


The open models include orientable and non-orientable space-forms which can 
be classified with glide reflections. On excluding glide reflection as a holonomy 
group, we get four orientable space-forms. The non-orientable space forms are 
not relevant to cosmology. The closed models on the other hand can be classified 
according to the different possible ways the opposite faces of a parallelepiped can 
be identified with each other. Another class of identifications can also be made 
on hexagonal fundamental polyhedron with rotations of 27/3 and 7/3. 


5.2 Spherical Manifolds 


Spherical manifolds have a universal covering of a compact hypersphere. A 3- 
sphere S? of radius R is the set of all points in 4-Dimensional Euclidean Space. 
The metric of the 3-sphere with coordinates x°, z!, z?, 3? can written as 


(x9)? + (z!)? + (z2)? + (23)? = R. (5.2) 


Converting to angular coordinates (x, 0, à), for x and 8 : [0, 1], 6 : [0, 27] 
x? = Rcosx, x! = Rsinxcos0, a? = RsinxsinÓcosó, x? = Rsin ysin ê sin ¢. 
We get the metric: 
do? = R?( dx? + sin x? (d6? + sin? 0d9?)] (5.3) 


from 
do? = (dx°)? + (dx)? + (dx?) + (da). 


The volume of the covering manifold is given by 


vol(S?) = f Av R? sin? yRdy = 21? P. (5.4) 
0 


Substituting r — sin x in metric, we get the FLRW metric form of the spherical 


manifold 
dr? 


dot = Pg 


+ r?(d6? + sinodo”) à; (5.5) 
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A good way of visualizing a 3-sphere is to use the analogy of a 2-sphere, 
where if we intersect the sphere with a plane and pass it through the sphere, the 
intersection points form circles of increasing diameter and subsequently of de- 
creasing diameter. Similarly, one can imagine the intersection of a 3-sphere with 
a 3-dimensional hyperspace and forming spheres of increasing diameter before 
reducing in size again to zero. 


The holonomy groups of S? were classified by [28] into subgroups T of SO(4) 
acting freely and discontinuously on S?: 


- Cyclic group of order p, Z,(p > 2) : Zp can be seen as generated by the 
rotations by an angle 27/p about some axis [0, ¢] of IR? 


- Dihedral group of order 2m, D,,(m > 2): Generated by rotations in the plane 
by an angle 2z/m and a reflection about a line through the origin. The operation 
preserves regular m-gons lying in the plane and centered on the origin. 


-Polyhedral Groups: Symmetry groups of the regular polyhedra in R? namely 
the Tetrahedral group of order 12, octahedral group of order 24 and Icosahedral 
group of order 60. The cube is included in the symmetry group of the octahe- 
dron and the dodecahedron is included in the symmetry group of the icosahedron. 


There is an infinite number of spaces that can be obtained by taking the 
quotent of S? with the above groups and varying the parameters p and m. The 
volume of the quotient manifold, M = S?/T obtained is given by 


Vol(M) = 22? R?/|T]. (5.6) 


5.3 Hyperbolic Manifolds 


Some of the most important contributions to locally hyperbolic spaces were made 
by Thurston [29], however these manifolds still remain much less understood than 
other homogeneous manifolds. Nevertheless, H? can be embedded in Minkowski 
space, IR? of metric 


ds? = —(dz')? + (da?)? + (dx)? + (dz)? 


as the hypersurface, 
-= (xt)? abs (27)? + (2°)? + t. 
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Thereby, the generators of the fundamental group G of H? can be related to ho- 
mogeneous Lorentz transformations [30]. 


Let us make coordinate transformations to introduce (x, 0, 9) with x € [0, inf), 0 € 
[0, 7], @ € [0, 27], we get 


x! = Reoshy, zr? = Rsinh x cos0, z? = Rsinh ysin cosó, z^ = Rsinh x sin@ sin ¢. 


We thereby get the metric for H? as: 
do? = R? ( dy? + sinh’ x(d0? + sin? 649?) . (5.7) 


This metric can be expressed in a more commonly used, FLRW form, of the 
metric, obtained by the coordinate change, r = sinh x that gives us: 


2 
do? = R? l T + dr? (d8? + sin? eas) . (5.8) 


The holonomies of H? can be described as the group of fractional linear trans- 
formations acting of the complex plane: 


ge O heated. 
cz+d 


This group is equivalent to the group of conformal transformations of R? which 
leave the upper half space invariant. 


The hyperbolic geometries for dimension > 2 are different from the 2-dimensional 
case. For instance, while a surface (genus > 2) can support an infinite number 
of non-equivalent hyperbolic metrics, a connected oriented manifold (n > 3) can 
only support at most one hyperbolic metric. This is confirmed further by the 
rigidity theorem which confirms that if two hyperbolic manifolds (dimension, 
n > 3) have isomorphic fundamental groups, they are necessarily isometric to 
each other. Hence for n > 3, the volume of a manifold and the lengths of its 
closed geodesics are topological invariants. 


For compact euclidean spaces, the fundamental polyhedron can possess only 
a maximum number of eight faces, despite allowance of arbitrary volume. In 
spherical manifolds, the volume needs to be finite and a fraction of the maximum 
possible volume S, i.e. S/T. For the case of hyperbolic manifolds however, there 
is no limit on the possible number of faces of the fundamental polyhedron. There 
is a lower limit however, on the minimum volume of the hyperbolic 3-manifold, 
a lower limit of which was set by Meyerhoff [31] to be Volmin > 0.00082 ?. 
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Chapter 6 


Evidence in favor of 
Multi-connectedness 


In 1900, Schwarzschild was already contemplating on whether the universe could 
have a finite and non-trivial topology. With the advent of modern technology, it 
has become feasible to test suggested theories and disregard or modify the ones 
that prove contradictory to observational evidence. The measurements of the 
temperature fluctuations in the CMB provided a solid observational basis to test 
current cosmological models. With the first set of data from the Planck satellite 
available, we have come a long way from COBE, which provided us with the first 
set of data quantifying perturbations in the microwave background. 


The observations from WMAP and Planck confirm the predictions made by 
COBE [32] as well as confirmed the observed anomalies in comparison with the 
expected results if the ACDM model is deemed correct. The anomalies were es- 
pecially evident on larger scales where the observed values of the multipoles were 
found to be much smaller than the theoretical values. This suppression of power 
has been widely explored in literature [34], [35], [36], [37], [38] as a possible evi- 
dence of the compactness of space. As pointed out in the introduction, a closed 
manifold would only allow certain harmonics to exist and the harmonics larger 
than the compact dimension will be suppressed, which would show in observa- 
tions as a suppression in large scale power. Therefore, it is possible that the CMB 
anomalies observed at large scales could be a signal of the non trivial topology of 
the universe. 


However the significance of these results as an invalidation of the concor- 
dance model is still under debate. It was shown by the WMAP team that the 
chances of the above anomaly occurring in a ACDM universe are one in 143 for 
the quadrupole and one in 666 for the overall large scale power. It has been 
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shown by Efstathiou [33] that these anomalies do not significantly deviate from 
the concordance model employing bayesian and frequentist methods as will be 
discussed below in section 6.3. Also, Oliveira et. al. [39] in their paper point 
out that the galactic cut that was used to remove foreground contamination from 
the map may have included some of the strongest hot and cold spots of these 
multipoles. 


However, a small universe is not the singular possible explanation for the 
anomalies. There are other models also which have ben put forward to explain 
the low large scale power, for instance a cutoff in the primordial power spectrum 
[41], possibly linked to the spatial curvature [40], multi field inflation models 
[42; 43] and quintessence models which involve a partial cancellation of the usual 
integrated Sachs-Wolfe effect [44]. 


6.1 Low Quadrupole 


The surprisingly small CMB quadrupole was first observed by COBE [45] and 
later confirmed by WMAP settling the initial doubts of the results being possibly 
noise induced. The results were further confirmed with the application of a re- 
fined analysis of Galactic foreground contamination. However, there still remains 
some uncertainty about the contribution of the cut in making the anomaly seem- 
ingly more striking. The observed quadrupole is the sum of cosmic and dynamic 
quadrupole and the contribution of the dynamic quadrupole can be determined 
from the dipole. A generic quadrupole has three orthogonal pairs of extrema (two 
maxima, two minima and one saddle point). However, the quadrupole that has 
been observed has its strongest pairs of lobes very close to the galactic plane, 
which could have removed a large fraction of the quadrupole power spectrum 
with the application of the galactic cut. The saddle point is close to zero and the 
quadrupole has a preferred axis in space along which the quadrupole has no power. 


The noise variance and beam issues are really low at large scales. Addition- 
ally the galactic foreground is only present in a small fraction of the sky and thus 
contributes negligibly to the total power spectrum. The WMAP team measured 
the quadrupole only from the part of the sky outside the galactic cut and hence 
the dominant uncertainty in its value was attributed to the foreground modeling. 
However, to get a better estimate of the effects of the galactic foreground, [46] 
measured the power spectra of the uncleaned ILC map and then subsequently 
zeroed out the power in the dirtiest regions successively leaving only the three 
cleanest regions. Band pass filtering of the resulting map produced almost the 
same spatial plots and power spectrum of the lowest poles as observed in the ILC 
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map, signifying that galactic modeling does not play a large role in the uncer- 
tainty of the multipoles. 


Using different combination maps and foreground cuts, we obtain slightly dif- 
ferent values of the quadrupole though they are all very low when compared to 
the expected theoretical value according to the ACDM model. A table compar- 
ing the theoretically expected and measured values obtained from independent 
methods is given below: 


Table 6.1: Comparision between the model and the measured values of the 
quadrupole and octopole [46]. 


Measurement óT2|uK^] p-value óT$[uK*| 
Spergel et al. model 869.7 855.6 
Hinshaw et al. cut sky 123.4 1.8% 611.8 
ILC map all sky 195.1 4.8% 1053.4 
Cleaned map all sky 201.6 5.196 866.1 
Dynamic quadrupole 3.6 


As evident from Table 6.1, Tegmark et. al. [46] obtained largely the same 
value for the quadrupole than the ILC maps. Thus even though Tegmark uses 
completely different modeling techniques, the quadrupole turns out to be virtually 
identical to the WMAP teams quadrupole which strengthens the argument that 
the contribution of galactic foregrounds to the quadrupole power spectrum is not 
significant. The p-value gives the probability that the given value of quadrupole 
is obtained given that the ACDM model is true and is calculated from a y? 
distribution for dT? with 5 degrees of freedom. 


6.2 Quadrupole-Octopole alignment and Planar 
Octopole 
Both the quadrupole and the octopole appear quite planar with their hot and cold 


spots centered on a single plane. To quantify the significance of the alignment, 
in [39], the angular momentum dispersion about an axis ù was defined by 


(| (AL)? |b) = SF maim (A), (6.1) 
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where aj, (fi) refers to the spherical harmonics in a rotated coordinate system with 
its z axis in the n direction, see Table 6.2. The axis around which this dispersion 
is maximized was then calculated. The preferred axis for the quadrupole and the 
octopole come out to be 


fa = (—0.1145, 0.5265, 0.8424), 


fi3 = (—0.2578, —0.4207, 0.8698), 
which is roughly in the direction of (l, b) ~ (—110°, 60°) towards the Virgo cluster. 


Table 6.2: Observed real values of quadrupole and octopole aj, coefficients in 
Galactic and Rotated(in direction of alignment) coodinates for the Cleaned map 
from [46]. 


m Galactic Rotated 
-2 -21.43 13.32 


l 

2 

2 -1l 6.03 0.40 

2 0 10.73 6.72 

2 1 -8.30 -0.40 
2 2 -19.39 28.86 
3 -3 40.71 50.58 
3 -2 2.45 -1.67 
3 -1 0.96 -0.50 
3 0 -6.52 -13.60 
3 1 -12.84 -0.27 
3 2 30.50 0.71 

3 2 -19.29 -20.68 


If the CMB is assumed to be an isotropic field, £9 and 73 should be indepen- 
dent of each other and subsequently the dot product between £9 and 3 should be 
randomly distributed between |-1,1]. In reality, the observed dot product between 
them comes out to be approximately equal to 0.9838. Thus, they are inclined 
to each other at an angle of 10.3?. The probability of an alignment this good 
happening is estimated to be 1 in 62 [39]. 


The octopole in the CMB has also been observed to be unusually planar, in- 


fact to a greater degree than the quadrupole. All other higher multipoles exhibit 
usual behavior as predicted by the ACDM model for a gaussian isotropic field. 
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The chances of the multipole being planar is 1 in 20 [39]. Oliveira et.al. [39] 
also defined a test statistic t, that measured the ratio of power from the m = |3| 
coefficients to the total power as 


"" mazs|a..a.a(&) |? + lass (ñ)? (6.2) 


3 x 
2 m3 asm (0)? 


This statistic for the rotated coordinates from Table 6.2 gives us a ratio of 94%. 
Monte Carlo simulations with independent values of aim coefficients corresponding 
to an isotropic Gaussian random field, we get higher t values than that calculated 
above only 7% of the time. Thus the contribution from |m| = 3 is dominant in 
the octopole. In comparison to the octopole, the quadrupole is not as significant 
in its planarity. 


6.3 Statistical Significance of the WMAP re- 
sults 


The significance of the low quadrupole and planar octopole are a topic of dis- 
cussion at present. Many factors are being taken into account which could have 
exaggerated the signal. For example, the effects of foreground contamination or 
alternative reconstructions of the true CMB signal from the five frequency band 
maps measured by the CMB might have contributed somewhat in skewing the 
original value. 


Inspite of the importance of paying attention to discrepancies observed in 
the observational test of an otherwise well confirmed model, we need to be sure 
that these discrepancies pass the significance tests of validity. Another ques- 
tion that may be asked is whether we get better overall fits to other modified 
models. These questions were explored by Efstathiou [33] in his paper where he 
applies the frequentist and bayesian analysis to arrive at a more realistic value 
for the chances of getting the same value, given that the ACDM Model holds true. 


It was shown by Spergel |47] that the observed values of the quadrupole and 
the octopole would occur with a probability of 1 in 700 if the concordance model 
was assumed to be correct. According to the current observations, there is an 
almost complete lack of signal at scales larger than 60?. Keeping this in mind, 
Spergel's [47] method was to calculate the statistic 
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Bes f " IC(6)f*d(cos), (6.3) 


for a large number of simulated CMB backgrounds for a ACDM universe. It 
was found that the probability of obtaining a value of S lower than the mea- 
sured value, was only 1.5 x 107%, thereby providing strong evidence that the 
observed values do not firmly support the ACDM universe. Efstathoiu [33], in 
his paper, used Monte Carlo Markov chains to calculate the posterior probability 
of the quadrupole given the CMB power spectrum and the covariance matrix 
and higher multipoles (800 < 1 < 2000) from CBI [50], ACBAR [49] and VSA 
[48]. He followed the MCMC analysis of [41] which used a six parameter ACDM 
model, namely a constant scalar spectral index n,, spectral amplitude A,, Hubble 
constant h, baryon density w = QA? , CDM density we = Qh? and redshift of 
reionization z,;; and evaluated the posterior distribution of the quadrupole for 
the given parameter set. 


The histogram of the resulting quadrupole amplitudes is shown in Fig. 6.1. 
The peak corresponds to the value, AT? = 1250 uK? with few values above 2000 
and below 1000 j,K?. This is quite a narrow distribution and relatively far form 
the observed value of AT? = 123K? in the publicly available data release. 


For a frequentist analysis, we take the distribution of C; which, for gaussian 
amplitudes of aim, will be given by a X? distribution, with further dependence on 
other factors such as the Galactic cut, instrumental noise and other sources of 
error 


RM | (21 1)81 dC; 
p m 


iP) = (Sr) (6.4) 


CF being the expectation value of Cj. Integrating Equation 6.4 gives us the 
probability of observing a value < C; as 
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Figure 6.1: A histogram of quadrupole values obtained for MCMC simulation 
universes for a ACDM model; the dashed line is the value of the quadrupole 
observed by WMAP [33]. 
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Efstathiou [33] applied the pseudo C; estimator to about 10° simulations with 
a galactic cut of 410? and obtained a histogram of the estimated quadrupole val- 
ues. Calculating the cumulative probabilities, the joint probability that AT? < 
123uK? and AT? < 61luK? was determined to be 0.3296 which is twice the 
probability calculated by Spergel et. al. [47]. With other small changes such 
as comparing the observed values with values from a lower allowed range of the 
theoretical spectrum or taking into account errors in the auto correlation func- 
tion that come from the inaccurate subtraction of the galactic foreground should 
decrease the significance of the observations even further. 


Efstathiou also applied bayesian analysis to the problem and calculated an- 
other set of significance level tests on the discrepancies between the observation 
and the model. A posterior probability analysis was done for the true values from 
the theoretical model manifesting a certain value on observation. According to 
the Bayes theorem, the posterior probability of a hypothesis H given data D is 
given by 


P(H/D) « P(D/H)P(H), (6.6) 


where P(H/D) is the probability of the hypothesis being true given a certain 
observation and P(H) is the prior probability of the hypothesis. Assuming the 
hypothesis that CT lies in the range CF + dC? , the posterior probability distri- 
bution for CT is given by 


1 (21+ 1)C; 
dP(CT —__._— parare d 6.7 
(€T) cerae (- Pa) acf (6.7) 
where C; is the observed amplitude and is also the maxima of the probability 
distribution curve of C7. 


Efstathiou [33] calculated the bayesian probability that the true values of 
the quadrupole and the octopole is greater that those predicted from the model, 
given the set of observations made (AT? = 123, K?, AT? = 611K? [51]). The 
discrepancies between the observed values and the peaks of the probability dis- 
tributions of the true values further reduce in this analysis with the ratios being 
p(C2)/p(CF) ria = 28 and p(Cs)/p(CZ)y;; = 1.6. Thereby, the frequency with 
which AT? > 1000uK?, given the observed WMAP quadrupole of 1234K? is 
only 0.10, which translates to only 1 in 10 chances. 
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Thus though the visible discrepancies suggest new physics, they do not entirely 
rule out the ACDM model. Keeping this in mind, we move on to discuss the 
topological models suggested to better incorporate current observations and the 
methods being devised to enable their detection. 


37 


Chapter 7 


Important Topologies 


Spatial curvature gives us a first guess at the class of topologies that could be pos- 
sible candidates to describe the universal manifold and can be directly calculated 
from an all sky analysis giving us a density parameter, ota Which measures the 
average mass energy density of space. As discussed in Chapter 2, depending on 
the whether the density parameter is less than, equal to or greater than 1, deter- 
mines whether the topology will be negatively curved, flat or positively curved 
respectively. Before dark matter was discovered, the total density parameter 
(Qio) only amounted to around 0.3, leading to a widespread and almost exclu- 
sive analysis of negatively curved topologies as possible models for the universe. 
At present, the value of the density parameter is estimated to be Quota, = 1.0240.2 
at the lo level [52] leading to a situation where any of the three curvature topolo- 
gies could be possible candidates, when taking observational errors into account. 


The current estimate was made by the WMAP team, whose data turned out 
to be well within the predictions for the concordance model spectacularly except 
on the largest scales (> 60°)[52]. The observed quadrupole is close to 0.15% of 
the expected value, while the octopole was less strikingly low but still about 7296 
of the expected value. A quadratic maximum likelihood analysis gives, however, 
somewhat larger probabilities in the range 3.2-12.5 % [33] 


If the universe were infinitely flat, the temperature fluctuations would have 
been expected to be present at all wavelengths just like a string of infinite length 
supports waves of all wavelengths. One of the arguments for the lack of such 
uniform presence of all temperature fluctuations is the compactness of space, 
thereby placing an upper limit on the wavelengths that can exist. Thus, cosmol- 
ogists now face the problem of finding a model for space that while producing 
the said anomalies on large scales, also closely follows the near perfect flat space 
predictions on smaller scales. 
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We study below some of the more widely discussed topological models that 
could produce the spectra that we observed in the CMB. 


7.1 Poincare Dodecahedron 


With the current estimates of the density parameter tilting ever so slightly to- 
wards a value larger than 1, a positively curved model seems to draw more favor 
as a viable topology. The simplest alternative was considered first, that of a 3- 
torus, which was indeed found to suppress the quadrupole. However, it was also 
shown to lower values of other multipoles, which did not mirror observations [53]. 
Some other models studied corresponded to lens spaces L(p,q), however these 
spaces suppressed the higher multiples even more than the lower multipoles, con- 
tradicting observations. 


It was observed soon that only well proportioned spaces selectively suppressed 
the lower multipoles. The spaces that fall in this category are the binary polyhe- 
dral spaces S?/T*, S?/O*, S?/I* which have positive curvatures and fundamental 
domains of a regular octahedron, truncated cube and a regular dodecahedron 
respectively. The Regular dodecahedron, commonly known as the Poincare do- 
decahedron was found to be within HRj,,./ R4, = 0.47 for Q = 1.02. The feature 
that set apart this particular topology from the numerous other positively curved 
models, that haven't already been ruled out due to obvious discrepancies, is that 
for a certain range of the density parameter, the poincare dodecahedron exhibits 
a strong suppression of power at l = 2 and a weak suppression at | = 3 when 
normalized to the l = 4 power. 


1,000 
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Figure 7.1: Comparision between multipoles from observation (light grey) and a 
flat (dark grey) and dodecahedral universe (black), from [54]. 
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The Poincare dodecahedral power spectrum depends on the mass energy den- 
sity parameter. It was found by Luminet et. al. [54] that the octopole power 
matches WMAP’s best predictions if (2,,, is between [1.010, 1.014], while for sup- 
pressed value of the quadrupole, Qot lies in the range [1.012, 1.014] as shown in 
Figure 7.2. Thus, it is quite clear that WMAP’s error range for the density pa- 
rameter comfortably includes both of these intervals, [1.02 + 0.2]. A comparison 
between the multipoles values obtained from the model to that from observation 
for Qm = 0.28 and Qa = 0.734 is shown in Figure 7.1. 


1.01 1.02 1.03 
ro 


Figure 7.2: Range of density parameter, (iota; for which the dodecahedral space 
agrees with WMAP observations, 1 and 3 being the observational data for 
quadrupole and octopole respective and 2 and 4 being the predictions from the 
dodecahedral model [54]. 


It is also striking that the Poincare dodecahedron does not have any degrees 
of freedom regarding its size or orientation. It is a positively curved space with a 
volume 120 times smaller than the simply connected hypersphere. It is fixed in 
its geometrical possibilities being a regular fixed volume dodecahedron and the 
only adjustable parameter is (,,4,. It is also globally homogeneous and thus looks 
the same to all observers at different positions. If we were to travel out of one 
side of the cell, we would return the same manifold albeit from the opposite face 
of the manifold. The adjoining edges (120?) of the spherical dodecahedron fit 
together snugly to tile the hypersphere. 


Yet, however appealing the striking correlations with the observation may 


seem, strict observational tests are a necessity for its consideration as a serious 
possibility, the matched circles test being a viable test. Cornish et. al. [55] pre- 
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Figure 7.3: A Poincare Dodecahedron tiled universe [2]. 


dicted detectable temperature correlations in the CMB for small multi-connected 
spaces. The radius of the six matched circles increases with increasing (3,4. If the 
density parameter were to be 1.013, our horizon radius being 0.38 in curvature 
radius units, then the in and out radii of the dodecahedron would be 0.31 and 
0.39 respectively, and hence the volume of the fundamental cell only 8396 of the 
observable sphere. The observable sphere in this case will intersect itself in six 
symmetrically placed circles of angular radius of about 35 degrees. However, no 
such circles were detected in the numerous independent searches made for the 
7-year WMAP data [57]. 


7.2 Bianchi universes 


One of the basic assumptions that the ACDM model relies on is that the uni- 
verse is homogeneous and isotropic on large scales. However, if we relax the 
assumptions on the global isotropy of space-time, we would get more generalized 
solutions of the Einstein’s field equations, giving us Bianchi universes. The re- 
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cently released topological analysis from the Planck team discussed Bianchi V Tj, 
models in detail, given that no statistically significant pairs of circles were de- 
tected in the searches performed on the high resolution data|11]. 


If we consider small anisotropy, we can get our solutions by applying small 
perturbations to the standard FRW model to evaluate the subdominant contribu- 
tion of the Bianchi component to the original temperature map from the ACDM 
Model [56]. 


The most general of the Bianchi Models are Bianchi VII, for flat, open models 
and Bianchi IX for closed models, and were analyzed in detail in [56] and [59] 
respectively. Bianchi VII, models were first compared to COBE in [62] and [63], 
and to WMAP in [60]. Jaffe et. al. [60] discovered a statistically significant 
correlation between one of the Bianchi VII, models and the WMAP ILC map by 
superimposing CMB from an unknown Bianchi component on the best-fit ACDM 
cosmology. Though promising, the parameters deduced from the Bianchi com- 
ponent did not agree with those from the ACDM. The important point to notice 
here though is that, on superimposing the Bianchi and ACDM data, some of the 
anomalies deduced in initial data disappeared. 


Later, a variation of the above analysis, including dark matter was explored in 
[61]. The results, as before still did not provide consistency between the param- 
eters derived from the Bianchi template and the CMB data. The first analysis 
of the scenario where the parameters of the Bianchi model are coupled to the 
parameters from the standard cosmological model was performed in [56]. 


Bianchi VII, models describe a universe rotating with an angular velocity 
w, and a three-dimensional rate of shear, specified by the antisymmetric tensor 
gij taken relative to the z axis. Throughout we assume equality of shear modes 
o = 012 = O43 [60]. Bianchi models contain a free parameter, x identified by 
Collins and Hawking[58] which further defines the parameter h for type VIL, 


models as [59] 
MEE 
x = 4| —.—. 7.1 
l= Qiotal l ) 


where the total energy density is the sum of matter density (Qm) and dark energy 
(Q4). x is related to the characteristic wavelength over which the principle axis of 
shear and rotation change orientation. It affects the tightness of the spiral-type 
temperature contribution to the CMB. Other defining parameters are the shear 
(c / H)o, vorticity (w/H) 9 and handedness « (-1 for left handedness and +1 right 
handedness). 
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The dimensionless shear (c;;/ H)o is thought to be related to the dimensionless 
vorticity, (w/ H)o as [59] : 


= EG, ea 


where the different shear modes have been assumed to be equal. The dimen- 
sional shear or vorticity characterize the amplitude of the induced temperature 
fluctuations but not its morphology. Apart from there being both left handed 
and right handed models, being anisotropic, the orientation of the Bianchi pat- 
tern is flexible on the sky. Thus it allows us three additional degrees of free- 
dom to define the model with, (a, 6,7) which becomes (0°, 0°, 0°) when the spi- 
ral pattern corresponding to Bianchi VII, model is centered on the south pole. 
Hence the Bianchi VII, models can be described by a seven parameter vector 
(Om, Qax, (w/ H)o, a. B, y) and & separates the left handed models from the right 
handed ones. 


The Planck team [11] simulated Bianchi CMB temperature maps using the 
BIANCHI2 code used by McEven et.al. [56] to simulate Bianchi VII, models. 
They simulated the Bianchi temperature maps for a varying x and density pa- 
rameters and the effect can be seen in Figure 7.4. 


For the most physically motivated open-coupled-Bianchi Model where the 
Bianchi VII, model is coupled to standard cosmology, there is no evidence in 
support of a Bianchi Contribution. The Planck team however does not completely 
rule out Bianchi VII, cosmologies in favor of ACDM cosmologies. 
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Figure 7.4: Bianchi universes withxe [0.1, 0.3, 0.7, 1.5,6.0] from top to bottom 
and (4, € [0.10, 0.30, 0.95] from left to right [11]. 
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Chapter 8 


Analytical Methods to detect 
Multi-connectedness 


A multi-connected universe consists of points which can be connected with more 
than one geodesic. Thus, an important observable effect of such a universe is 
that light will have travelled across the fundamental domain more than once if 
its size is lesser than the surface of last scattering. For a small enough universe, 
we should be able to see multiple images of a cosmic object in different directions 
in the sky. Looking for these multiple images is one of the most straightforward 
and direct methods of looking for evidence of multi-connected topology. 


The possibility to beyond the fundamental cell of a multi-connected universe 
depends on the odds of it being comparable to the size of the observable universe. 
The chances are maximized when the search is conducted on the farthest possible 
surface that we have access to data from. The cosmic microwave background offers 
the best possible chances with current methods to detect a topological signature 
which would skew the otherwise isotropic nature of the temperature perturba- 
tions. In other words, the temperature perturbations in direction £4, T(f£i) would 
be correlated to the perturbations in another direction m by an amount that 
would clash with the one expected from a usual isotropic correlation function 
C(@) defined in section 8.1, 0 being the angle between ñ and rh. In a pixelized 
map, this results in a correlation index than does not solely depend on the angular 
separation between pixels. Such correlations are an important clue to detecting 
possible manifestations of topology. 


If we live in a multi-connected universe, with the fundamental cell size smaller 
than the observable size of the universe, we should in principle be able to detect 
multiple images of the most distant cosmic objects. However, a direct search for 
these multiple images brings us face to face with many observational problems 
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like: 


(i) Morphological effects when viewing from different directions, 


(ii) Differentiating or finding correlations between objects at different dis- 
tances which would correspond to different stages in the life of the object and 


(iii) Foreground contaminations or high obscuration regions masking the ob- 
ject or its images. 


These problems, however, can be overcome by using statistical methods, which 
provide a more accurate measure of the possible matches by using indicators or 
signatures to search for a given topology. 


Another way to look for non-trivial topology is to use data from the cosmic 
microwave background which gives us information about the pressure density 
fluctuation map of the universe at the surface of last scatter. Many methods 
have been suggested to look for topology using the CMB and can be looked up in 
detail from [66]. The most popular from among these is a method that looks for 
matching circles in the background, called the circles in the sky method, devised 
by Cornish and Spergel in 1998 [55] and is further discussed in Section 8.3 


8.1 Correlation Matrices 


CMB temperature correlation matrix is a double radial integral of the average of 
the product of source functions that represent the transfer of the photons through 
the universe from the last scattering surface to the observer: 


Xlss Xlss 
y= f dy f dy! (S(x4p) S dy). (8.1) 
0 0 


where q, and d, are unit vectors in the direction of the pixels p and p’ while x 
and x’ are proper distances along radial rays towards the last scattering surface. 


There are currently two ways of computing CMB correlation functions for 
multi-connected universes. One of the ways to compute the CMB correlation 
functions for multi-connected universes would be to construct an orthonormal 
set of basis functions that satisfy the boundary conditions of the fundamental 
domain of the manifold. One then computes the spatial correlation function 
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(S(xG>)5(x'Gp')) from this basis [55; 64]. 


The usefulness of these methods relies on the fact that in case of compacti- 
fication and/or global anisotropy, the correlation function would not solely be a 
function of the angular separation between p and p’ . In harmonic space, the two 
point correlation function is given by 


mm' — (i Oe) Z Clómm OW ; (8.2) 


where am are the spherical harmonic coefficients of the temperature function on 
the CMB given by 
T(G) = V ai Yos). (8.3) 


im 


The Planck collaboration [11] calculated these correlation functions for differ- 
ent topologies and plotted them to obtain the correlation patterns typical of a 
certain topology. Figure 8.1 shows the correlation matrices in pixel space, show- 
ing the magnitude of correlation between the pixels. Smoothing with a gaussian 
field gives us a microwave background map for each of the said topologies as 
shown in Figure 8.2 


8.2 Cosmic Crystallography 


The most evident way to look for signatures of a multi-connected topology is to 
look for repeated images of distant cosmic sources. This search however is not 
without difficulty and involves problems such as having to take into account the 
time evolution of the object which may look considerably different from its image 
at different distances in the cosmic scenario. Morphological and foreground dis- 
tortion could further skew the data and make identification of images harder that 
it would initially appear. It seems best to approach the problem with a statistical 
approach by using specific statistical signatures to look for signs of a non-trivial 
topology. Here, we describe the pair separation method, more commonly known 
as the Cosmic Crystallography method as described in detail by [26]. 


In 1996, Lehoucq et. al. [65] proposed the Pair Separation Histogram (PSH) 


method to look for correlations in the sky. The method relies upon the func- 
tion of the distance, F(d) = s between pairs of images in a catalogue C. Then 
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Figure 8.1: A correlation structure in pixel space of i) At the top, a simply 
connected isotropic universe ii) second row, from left, dodecahedral, octahedral 
and equal sided torii, with the fundamental domain comparable to the size of the 
last scattering surface iii) Subsequent rows correspond to a smaller fundamental 
domain with respect to the last scattering surface, from [11]. 


Figure 8.2: Realizations of the microwave background from the correlation maps 
superimposed with a gaussian profile of full width at half maximum of 640’, from 
[11]. 


we calculate the number of pairs whose separation lie within one of the sub in- 
tervals J;, from the partitions of [0, Smar] where Smar = F(2Xmaz) and Xmas 
is the survey depth of the Catalogue. The separation function is usually taken 
to be simply the distance between the pair s = d or its square s = d?. A nor- 
malized plot of the separations gives us the pair separation histogram of the data. 


In more detail, for a catalogue of n cosmic sources, 7(s) denotes the number 


of pairs of sources whose separation is s. Divide the interval (0, Smar] in m equal 
subintervals(bins) of length ôs = Smar/mM; the intervals are thus 


—L.i-21,2,..,m (8.4) 


and centered at s; = (i — 1)ós. The Pair Separation Histogram (PSH) will be 
given by the counting function 


(8.5) 


Figure 8.3: A PSH for a 3-torus topology where the horizontal axis is the sepa- 
ration squared while the vertical axis plots the number of pairs [26]. 


An example of a PSH for a universe with non-trivial topology, a 3-torus in 
this case, is shown in Figure 8.3. The curve characteristics are a mean curve with 
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additional spikes which are a combination of noise and the signal. The mean 
curve is representative of the expected PSH, EPSH and is exclusive of statistical 
noise. 


Further tests revealed confirmation of expected spikes in flat manifolds [65], 
yet no spikes were observed for hyperbolic multi-connected manifolds [67]. The 
reason was discovered after a theoretical statistical analysis, that the spikes of 
topological origin in PSH were due to only one type of isometry, the Clifford 
translations (CTs) which are isometries g € I such that for all p € M, the dis- 
tance d(p, g:p) is a constant. The CTs reduce to regular translations in Euclidean 
spaces and the lack of CT translations in Hyperbolic geometry results in its ab- 
sence of spikes in these universes. 


Though other isometries also result in a slight deformation of the expected 
pair separation histogram as shown by [68], their amplitude is usually small as 
compared to the background noise and thus, they are barely detectable. 


Although a good method for detecting non-trivial topology for flat mani- 
folds and a few spherical manifolds that admit Clifford transitions, the above 
method fails for hyperbolic universes. Uzan et. al. [69] thereby suggested an- 
other method called the Collecting Correlated pairs method which uses the 
property of preservation of distances between images. 


The method's central premise is that the distance function between two ob- 
jects and their images should be equivalent, i.e. for each g € T, 


d(p, q) = d(gp, gq) (8.6) 


given that the pair of images(p,q) is also in the Catalogue. The analysis is 
basically done computationally by arranging the distances between all P — n(n— 
1)/2 possible pairs between the n sources in the catalogue and arrange them in 
increasing order. Then, the difference between successive distances is evaluated, 
with N the number of times the distance equals zero or a very small number e 
taking measurement errors into account. The CCP index, which is given by 


N 
R = —, 8.7 
Pol (8.7) 
is evaluated. A value of R > 0 can be taken as a sign of multi-connectedness 
with higher values making the case more certain. The method however, has its 
own associated problems which comes from uncertainty in the determination of 
precise placements and distances of cosmic objects, not to forget the red shift 
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corrections due to additional peculiar velocities of these objects. [70] discusses 
further limitations and uncertainties associated with the above methods. Nu- 
merous other methods are currently being developed which may provide a more 
accurate statistical indicator for non-triviality [71; 72; 73]. 


8.3 Matched Circles Method 


The Cosmic Microwave Background is currently the most primitive source of in- 
formation that we have access to and it allows us to look as far into space as 
possible with current technology. The information in the CMB encodes the data 
from redshift, z = 1100 which is the last scattering surface, the radius of which 
is given by 


2-0 
Reis = Rar, arccosh ( J ; (8.8) 
Qo 


where the photons first decoupled from matter. 


The main premise of this method relies on the idea that if the observable 
universe around us can be thought of as a sphere of radius xy;,4. In the scenario 
where the fundamental domain of our topology is smaller than this sphere, there 
will be an intersection between the spheres of their last scattering surfaces, and 
the intersection would be the shape of a circle irrespective of the curvature or 
structure of the topology. However, the positions and angular correlations of the 
circles in the sky will be determined by the topology. The matched circles will 
have an identical patter of temperature variations. The mappings from the sur- 
face of last scatter to the night sky is conformal thus making sure that the angles 
are still preserved. This implies that the identified circled appear equally resealed. 


Cornish et. al. [55] developed the statistical method to detect these correlated 
circles in the microwave background. They selected two points, p and r in the 
night sky and drew circles of angular radius a around each point and consider all 
possible phases, ¢, between the points. If T,(¢) and T,.(¢@+¢.,) are defined as the 
temperature fluctuations along the two circles, the circle comparison statistic is 
given by 
Í (T,(6)? + T(¢)?) 


where () = DA dó and the + sign refers to whether the circles are ordered in the 
same or opposite directions respectively. S varies from [-1,1], being + for perfect 
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correlation and 0 for perfectly unmatched circles. For each radius a, we com- 
pute the maximum possible value of the S statistic. The circles are ordered in 
clockwise direction for orientable and anti-clockwise for non-orientable topologies. 


Better angular resolution provides a better possibility to detect these corre- 
lations as it provides a larger dataset of points. An experiment with angular 
resolution A0 provides around N ~ 27sina/Aé@ data points around each circle 
of angular radius a. 


Cornish et. al. [74] later modified the statistic such that the above compu- 
tation can be realized faster by working with Fast Fourier Transforms (FFT) of 
the temperature correlations along the circles, 


Tol) = M ^ Ty exp(imó) 


where |m| which represents the mth harmonic around the pth circle. The statistic 
can be written as 


JAG $*) = 2c Smerp(—imds) 


where 


Sm = 2 » al al yuQ A + Mal) 


using the inverse fourier transform. The modified statistic can be written as: 


25n ca re eg Pes 


+ — 
9 (e$) = Dn a 


(8.10) 


which weighs the temperature coefficients with the mth harmonic taking into 
account the number of degree of freedom per mode. Doing so, enhances the con- 
tribution from the small scale structure as compared to the larger fluctuations 
which are dominated by the Integrated Sachs Wolf effect (ISW). 


General searches explore a six parameter space: the location of the first circle 
centre (0p, ¢,), the location of the second circle centre (6,,¢,), the angular radius 
of the circle a and the relative phase between the two circles ¢,. Being excessively 
computationally intensive, with N? logN operations, the searches performed till 
now, have been limited to antipodal or nearly antipodal circles. There is a large 
subset of topologies that could manifest in an antipodal configuration, except 
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perhaps the Hantzsche-Wendt space for flat universe or the Picard space [75]. 


The searches performed by the Planck team turn up no evidence for a multi- 
connected topology [11]. There were no statisticaly significant correlations above 
the false detection level of Amin ~ 20°. Thus any topology predicting matching 
pairs of back to back circles larger that amin have been ruled out to a confidence 
level of 99%. 


8.4 Bayesian Analysis 


Though observational methods provide a concrete basis for ruling out or confirm- 
ing a certain topological model, they are rife with measurement errors. It is thus 
useful to study the posterior likelihood of various models theoretically before- 
hand. Especially with the high resolution data from Planck released, searches 
for matched circles by the Planck team [11] did not provide any statistically 
significant match to any topology, thus making introduction of other methods 
mandatary. 


These methods take into account that Gaussian likelihood functions apply as 
much to anisotropic and topological models and to the standard cosmology, albeit 
without isotropic correlations. The form of the likelihood function depends on 
the anisotropic model. For instance, for multi-connected models, the pixel-pixel 
correlation matrix includes anisotropic correlations while Bianchi models give a 
deterministic template in addition to the standard isotropic correlations for stan- 
dard cosmology. It is assumed that there are no other sources of non gaussian 
signals. 


Let us denote the model under examination by M, which could be an isotropic 
universe model, a Bianchi model, or a non-trivial topology. The data is denoted 
by vector d, containing harmonic coefficients or pixel temperatures depending 
on the basis functions used for representing the temperature function. Vector 6 
are the cosmological parameters which define M, and they can be split into two 
sets, 9¢ which are parameters shared with the isotropic and simply connected 
model and 04 which are the parameters specific to the anisotropic properties of 
the model. The posterior probability estimation of the parameters uses the Bayes 
theorem as given by 


P(8/M)P(d/0, M) 
P(dM) ^ 


P(0/d, M) — (8.11) 
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where P(@/M) is the joint prior probability of 04 and 0c, P(d,0, M) = £ is the 
likelihood. 


The method is quite similar to the methods used in standard cosmological pa- 
rameter estimation except that the above method is much more computationally 
expensive, given the additional complexity. 


For the analysis of non-trivial topologies, the 0c are the parameters that are 
shared with the concordance model while 04 are the parameters specific to that 
particular topology for instance compactification lengths, the curvature parame- 
ter etc. Thus, it is usually 0r parameters that have to be varied for different analy- 
sis. The 2-point signal correlation matrix, Cpp is precomputed for the topologies. 
Parameters such that signal amplitude, A and the orientation of the fundamental 
domain with respect to the sky, o are usually maximized or marginalized. 


Having defined the parameters, we define the Likelihood, which the probability 
that we would get a certain data vector, d along with its noise matrix, N, given 
a certain non-trivial topology. It is expressed as [11] 


P(d|Cl6c, 6r, T], A, d) « exp[- I (AC +N) dg}. —— (812) 


1 
AC +N 


Usually, likelihood calculated in pixel space is not entirely accurate and all 
modes are not equally significant for the analysis. Thus, the quantities dp, Cop, Npp 
are expanded into discrete mode functions w,,(p) orthonormalized over the pix- 
elized sphere, tentatively with weights w, such that 


X w(pyus (p) (p) = 9v 


'The role of the weights is in improving the accuracy of transforms on a pixelated 
sky. Now, the coefficients of the quantities in question can be written as: 


d, = > dy; (p)w(p), (8.13) 
Con = SOS) Coyne (v w(p)w(p), (8.14) 
Naw = >> >> Cop talpi (p )w(p)w(p). (8.15) 


54 


Nm modes are then selected out the series for comparison with the real signal, 
while other modes are marginalized. The likelihood equation now looks like this: 


N 
1 m 
ezpi-; J d (AC + N).Ld,). (8.16) 


1 
V|AC t N|u nl 


p(d|C 0c. 0r. T), V, A) x 


The quantities C and N above only include the chosen modes Nm. The choice 
of the mode functions and their number is a trade off between the ease of invert- 
ing the C+N matrix and the useful information included in the data. 


In the case of determining likelihood for Bianchi cosmologies, the analysis 
involves fitting a certain deterministic Bianchi template on a stochastic CMB 
background [56]. The parameters are denoted by 0c for the stochastic back- 
ground with a power spectrum Cj;(0c) and 0g for the Bianchi template. The 
likelihood of the model is denoted by 


1 2 
P(d|85, 06) ex JE * (9c, 0B)/2], (8.17) 
where 
X (0c, 0B) = [d — t(65)] X ^ (6c)|d — &(65)] 


and d = {dim}, b(05) = {bim(@s)} are the spherical harmonic coefficients of the 
data and the template respectively, evaluated unto Ima, = 32 which was deter- 
mined to be enough to capture the structure of the CMB [76]. 


Incorporating the rotations to the Bianchi template, the likelihood is evaluated 
inharmonic space. The Covariance matrix X(06) is determined by the mask 
settings and is equal to C(0c) without any mask applied, where C(0c) is the 
diagonal CMB covariance matrix [77]. With a partial sky mask applied however, 
we have 

X(0c) = C(0c) + M 


where M is a non diagonal covariance matrix for the mask. When the Bianchi 
template is null, the likelihood equation 8.17 reduces to the expression used for 
usual cosmological parameter analysis for the concordance model. 


The Planck team [11] used two separate ways to do the likelihood analysis for 


the template, in terms of the coupling between the Bianchi template and stochas- 
tic background parameters 05 and 0c. Though the coupled case is favored over 
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the decoupled case by nature, it is worthwhile to examine the other case as well. 
For the decoupled case, they took a flat cosmological model while for an open 
universe was used in the coupled case for consistency with the Bianchi VIT, mod- 
els. 


For each of the models, one can evaluate the Bayes factor to examine their 
credibility over the concordance model, from 


E = P(d/M) = f dóP(d/0, M)P(0/M). (8.18) 


If the prior probability of the model is unknown, the Bayes factor can be 
determined by the ratio of the model’s likelihoods. Jefferys [78] proposed the 
Jeffery’s scale which employs the log Bayes factor to compare different model 
likelihoods. It is given by 

A ln E = In( E/E) 


and is a measure of the relative likelihood of models E, and E». A log Bayes 
factor AlnE > 5 is regarded as conclusive evidence of the superior significance of 
model E, over E; (given that E, > Es). The certainity of the same reduces for 
a smaller Bayes factor with A In E < 1 considered inconclusive. 
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Chapter 9 


Conclusions 


The primary motive of this thesis has been at providing a comprehensive review 
of the current focus of Cosmic Topology. Though the underlining principles have 
remained largely unchanged through the active years in the search for a defining 
topology of the universe, the focus of the searches has varied significantly. A 
crucial beginning point in topology is the discussion of the fact that Einstein’s 
equations of general relativity only constrain local properties of space-time and 
do not specify in any way global properties of the universe. 


One of the most evident changes from my perspective was the shift in focus 
from hyperbolic topologies during the time of COBE (before the discovery of 
dark energy), which provided evidence (Oyo: = 0.3) that the universe might have 
a negative curvature, to now, where the focus is on positively curved topologies, 
the Poincare Dodecahedron being the most notable, when the odds lie more to- 
wards a positively curved topology (Qio = 1.02), though not entirely ruling out 
a hyperbolic one. Another noticeable evolution has been in the kind of methods 
that have been used to investigate multi-connected topology. Having started out 
as a purely theoretical speculation, cosmologists gradually devised and developed 
observational tests as the skies were mapped to an unprecedented degree of res- 
olution and it became possible to measure spatial curvature, spectral index etc 
thereby quantifying parameters for quantitative analysis of various topologies. 
This enabled cosmologists to put reliable constraints on cosmological parameters 
for different topologies and test them against observations. 


The theoretical methods focussed on by the Planck team were used to quan- 
tify relative likelihood of topologies and garner theoretical evidence in favor of 
or ruling out different topologies. Identifying in advance the topologies which 
agree to a reasonable degree with parameters obtained from observation, gives 
us a better chance at detecting those topologies by looking for specific signatures 
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unique to the topology. The methods used by the Planck team, be it evaluating 
likelihoods or searching for circles in the sky, did not succeed in providing any 
conclusive evidence for any particular non trivial topology. 


The universal manifold is often likened to a musical instrument, which only 
plays certain harmonics thereby offering a possibility to know its shape and size 
from merely observing the sounds it produces. The strings here are analogous 
to the photon baryon plasma at the surface of last scatter and the sound to the 
acoustic oscillations present in the plasma. A straightforward inference of this is 
that the harmonics of the oscillations that existed in it had to be limited by the 
shape of the manifold they existed in. That these harmonics should still be de- 
tectable from the CMB photons that emanated from the plasma, is a reasonable 
assumption, given that it is not masked by foreground or perturbations collected 
on the way. 


The possibility of the anomalies of the low quadrupole and planar octopole 
being manifestations of a non-trivial topology is a promising evidence of the 
premise that there may be detectable local signatures of the universal manifold. 
However until some more ingenious methods, like the brilliant circles in the sky 
method, are developed to determine our topological envelope measurable from 
earth or somewhere in its proximity, we can only rely on technological enterprise 
to further our investigations. Future missions to map the microwave background 
should be able to provide us with higher resolution data with a better signal- 
to-noise ratio which will further rule out some topologies and bring others into 
perspective. 
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